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ABSTRACT 

Predictions of turbine vane and endwall heat trans- 
fer and pressure distributions are compared with exper- 
imental measurements for two vane geometries. The 
differences in geometries were due to differences in the 
hub profile, and both geometries were derived from the 
design of a high rim speed turbine(HRST). The ex- 
periments were conducted in the Isentropic Light Pis- 
ton Facility(ILPF) at Pyestock at a Reynolds No. of 
5.3 x 10 6 , a Mach No. of 1.2, and a wall-to-gas tem- 
perature ratio of 0.66. Predictions are given for two 
different steady state three-dimensional Navier-Stokes 
computational analyses. C-type meshes were used, and 
algebraic models were employed to calculate the turbu- 
lent eddy viscosity. The effects of different turbulence 
modeling assumptions on the predicted results are ex- 
amined. Comparisons are also given between predicted 
and measured total pressure distributions behind the 
vane. The combination of realistic engine geometries 
and flow conditions proved to be quite demanding in 
terms of the convergence of the CFD solutions. An ap- 
propriate method of grid generation, which resulted in 
consistently converged CFD solutions, was identified. 

N omenclat ure 

c - True chord 

c r - Axial chord 

d - Distance from surface 

k - Thermal conductivity 

Mis - Isentropic Mach No. 

Nu - NusseJt No. based on true chord and k(T{) 

P* - Total pressure 

Pr - Prandlt No. 


Pr t 

- 

Turbulent Prandlt No. 

Re 2 

- 

Exit true chord Reynolds No. 

s 

- 

Surface distance 

V 

- 

Total temperature 

Tu 

- 

Turbulence intensity 

U 

- 

Velocity 

vt 

- 

Normalized distance of first grid line 

e 

- 

Momentum thickness 

vt 

- 

Turbulent eddy viscosity 

p 

- 

Density 

Subscripts 



EXIT 

- 

Exit of computational domain 

INLET 

- 

Inlet of computational domain 

M 

- 

Measurement plane 

S 

- 

Surface of blade 

1 

- 

inlet, or surface 

2 

- 

outlet, or surface 


INTRODUCTION 

Confidence in the validity of three-dimensional 
Navier-Stokes predictions for turbine aerodynamic per- 
formance and surface heat transfer characteristics in- 
creases when computations are shown to predict mea- 
sured values for cases which approach actual engine ge- 
ometry and flow conditions. For example, the engine 
Reynolds and Mach numbers, as well as the turbulence 
intensity should be matched as closely as possible. Con- 
fidence levels increase when the test geometry incor- 
porates features which axe part of a three-dimensional 
design philosophy. A turbine vane design philosophy 
could include convergence in the meridional plane, as 
well as underturning of the vane in the endwall region 
to reduce secondary losses. When design considerations 
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such as these are present, blade and endwall heat trans- 
fer may be substantially different from values measured 
in either a linear or annular cascade with constant sec- 
tion vane geometry. 

Comparisons with experimental turbine blade heat 
transfer data which include both endwall and blade 
results provide a significant test of three-dimensional 
Navier-Stokes analyses, since the endwall heat trans- 
fer is significantly affected by secondary flow patterns. 
The blade heat transfer is more two dimensional in na- 
ture. Blair(1974), and Boyle and Russell(1990) pre- 
sented endwall heat transfer distributions obtained in 
linear cascades for large scale vanes, at relatively low 
Mach numbers. York et ad. (1984) obtained similar data, 
but at transonic Mach numbers. Arts and Heider(1994) 
presented vane and endwall heat transfer measurements 
obtained at transonic Mach numbers in an annular cas- 
cade with constant radii, endwalls. 

Test cases with hubs and shrouds of constant radii, 
and with blades of little or no twist provide useful in- 
formation for predicting aerodynamic and heat transfer 
behavior. However, turbine blading employed in actual 
engines typically has vanes with significant endwall con- 
vergence and twisted blades. Comparing predictions 
for test cases with these geometric features increases 
confidence in the analysis’ ability to predict aerody- 
namic and heat transfer performance at actual engine 
conditions. The previously cited experimental data are 
for geometries with constant radii hub and shroud sur- 
faces, and with little or no blade twist. A series of ex- 
periments measuring turbine aerodynamics and surface 
heat transfer utilizing a vane geometry representative 
of an actual engine have been conducted in the DRA 
Isentropic Light Piston Facility(ILPF). Kingcombe, Ha- 
rasagama, Leversuch, and Wedlake(1989), presented 
the overall performance and detailed vane heat trans- 
fer and aerodynamic measurements for a transonic tur- 
bine. Vane measurements were given for a range of 
Reynolds numbers at an exit Mach number of 1.14. Ha- 
rasagama and Wedlake(1991) presented experimental 
surface pressure and heat transfer distributions on the 
vane and on the hub and casing endwalls for a range 
of Reynolds and Mach numbers. Both of these ref- 
erences were for a vane with reduced turning at the 
hub and tip casing, but with constant radius endwalls. 
Chana(1992), measured vane, hub, and tip heat transfer 
and pressure distributions using the same vane geome- 
try for two different hub contour geometries. 

In recent years a number of three-dimensional 
Navier-Stokes predictions for vane and/or rotor heat 
transfer have been presented. Choi and Knight(1988), 
Hah(1989), and Ameri and Arnone(1994) showed com- 


parisons with large scale-low speed heat transfer data 
of Graziani et al.(1980) using both algebraic and two 
equation turbulence models. Chima et al.(1993) devel- 
oped an algebraic turbulence model, which was used 
to compare vane endwall heat transfer predictions over 
a range of Reynolds numbers with the experimental 
data of Boyle and Russell(1990). In addition to predic- 
tions for isolated blade rows, three-dimensional Navier- 
Stokes heat transfer analyses have been performed for 
entire turbine stages. Boyle and Giel(1992), and Ameri 
and Amone(1992) showed comparisons with the exper- 
imental heat transfer data of Dunn et al.(1994) and 
Blair(1994). Heider and Arts(1994) showed compar- 
isons of predicted heat transfer with the data of Arts 
and Heider(1994) using the Baldwin- Lorn ax turbulence 
model. Except for the comparisons of Heider and 
Arts(1994), these comparisons have been for subsonic 
conditions. It remains to be seen if algebraic models 
can predict turbine heat transfer as accurately as two 
equation models for high Mach and Reynolds number 
cases. 

The ideal turbulence model would give accurate 
heat transfer predictions with little computational over- 
head. Algebraic turbulence models require less com- 
putational effort compared to multiple equation mod- 
els. They also have the advantage in that empirical 
correlations can be easily incorporated into the model- 
ing to account for effects such as transition and aug- 
mentation of the leading edge heat transfer due to 
freestream turbulence. However, if the algebraic mod- 
els are not accurate, their use is not warranted. Ameri 
and Arnone(1994) showed heat transfer predictions us- 
ing algebraic and two-equation models. Their results 
indicated that the degree of agreement with the ex- 
perimental data was not significantly improved using 
a two-equation turbulence model. Results using either 
turbulence model were in good agreement with the ex- 
perimental data. Other factors, such as the thermal 
boundary condition and spanwise grid density, resulted 
in heat transfer variations as large as between the two 
turbulence models. Algebraic models also require less 
CPU time per time step, and less core memory than 
higher order models. Ameri and Arnone(1994), showed 
that heat transfer predictions using two-equation mod- 
els required nearly twice the CPU time to converge, 
compared with an algebraic model solution. The high 
Reynolds number cases investigated herein require large 
size grids, which favors algebraic models. Also, use 
of higher order models would have required additional 
computer simulations to determine the effects of the as- 
sumed value of the inlet length scale on the predicted 
heat transfer. The length scale at the inlet was not 
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known, and Ameri and Arnone(1992) showed that blade 
surface heat transfer is significantly affected by the as- 
sumed value for the length scale at the inlet boundary. 

As the result of a cooperative agreement between 
NASA Lewis Research Center and the Defence Research 
Agency(DRA, Pyestock) a study was undertaken to 
predict the turbine vane heat transfer of two geomet- 
ric configurations. The predictions are compared with 
the experimental measurements of Chana(1992). These 
measurements included pressure and heat transfer dis- 
tributions on the vane surfaces as well as the hub and tip 
casings. Pressure and heat transfer distributions were 
measured on the vane, hub, and shroud surfaces for two 
different hub endwall geometries. In addition, the total 
pressure distribution behind the vane was measured for 
one configuration. These data are for a high Reynolds 
number of Re 2 = 5.3XI0 6 , and at a pressure ratio cor- 
responding to an exit Mach number of 1.2. The high 
Reynolds number requires small grid spacing in the near 
wall region. This, in turn, results in large computa- 
tional grids. Because of the large computational grids, 
and indications that algebraic turbulence models would 
give accurate heat transfer predictions, algebraic turbu- 
lence models were used. Heat transfer comparisons are 
shown using three different models for turbulent eddy 
viscosity. Predictions were made for the vane surface, 
hub, and casing heat transfer. Also, the ability of the 
analyses to predict surface static pressures, and the to- 
tal pressure distribution behind the vane is examined. 
The reason for comparing predicted wake profiles with 
measurements is that these comparisons may shed light 
on the differences in the heat transfer results. Good 
agreement in the downstream pressure distributions in- 
dicates that the analysis predicted secondary flows cor- 
rectly. 

DESCRIPTION OF COMPUTATIONAL 
ANALYSIS 

Two different steady state three-dimensional 
Navier-Stokes computer code analyses were used. One 
computer code, referred to as TRAF3D, is the finite 
volume analysis described by Arnone et al.(1992), The 
other code, referred to as RVC3D, is the finite difference 
analysis described by Chima(I99I) and by Chima and 
Yokota(1988). The reason for using two Navier-Stokes 
analyses is to help insure that the conclusions drawn 
from this work are applicable to a range of CFD analy- 
ses, and that they are not specific to any one formula- 
tion. While the two Navier-Stokes analyses had differ- 
ent approaches to the discretization, they were similar 
in other respects. Both analyses used a Runge-Kutta 


scheme to march the solution in time, and both em- 
ployed implicit residual smoothing. 

The code TRAF3D uses the turbulence model de- 
veloped by Baldwin and Lomax(1978). The RVC3D 
code uses either of two turbulence models. One is a 
variation of the Baldwin-Lomax model, and results ob- 
tained using this model are labeled as Chima’s model. 
The other is a variation of the Cebeci-Smith(1974) 
model, and results obtained using this model are la- 
beled as Cebeci-Smith results. Chima et al.(I993) dis- 
cuss the implementation of these two models. Unless 
stated otherwise references to the TRAF3D code im- 
ply the Baldwin-Lomax turbulence model, and refer- 
ences to the RVC3D code imply Chima’s turbulence 
model. Both of these models use the Baldwin-Lomax 
transition criterion. The transition criterion in the 
Baldwin-Lomax model does not account for the effects 
of freestream turbulence. Since these tests, and ac- 
tual engine operation, are at relatively high levels of 
turbulence, the transition criterion of Mayle(1991) was 
incorporated into the turbulence model. Also, high 
freestream turbulence results in leading edge Frossling 
numbers significantly greater than unity. The model of 
Smith and Kuethe(1966) was used to account for the 
effects of freestream turbulence. Pr and Pr t were held 
constant at 0.70 and 0.90 respectively. 

The inlet total pressure was 5.0 bars. The inlet 
total temperature was 440° A', and the surface temper- 
ature was 290° K. The hub static- to-total pressure ratio 
was 0.42. The Reynolds number based on exit condi- 
tions, and true chord was 5.3 x 10 6 . A uniform tem- 
perature boundary condition was imposed on all solid 
boundaries. Span wise radial equilibrium was assumed 
at the exit boundary of the calculation domain. At each 
spanwise location the exit static pressure was allowed 
to vary in the circumferential direction. The average 
static pressure was specified at the hub exit, but the 
pitchwise variation was determined from the internal 
flow field. The hub and tip inlet full boundary layer 
thicknesses were determined from measurements to be 
9% of axial chord. Uniform toted conditions were spec- 
ified for the inlet core flow, and uniform static pres- 
sures were specified through the inlet boundary layers. 
The inlet temperature and velocity profiles through the 
boundary layers were determined using flat plate cor- 
relations. The friction factor and Stanton number are 
first determined from the specified inlet boundary layer 
thickness. Once these quantities are found, the profiles 
are determined from correlations given by Kays and 
Crawford(1980). Using a simple power law for the inlet 
temperature profile produces an erroneous result, since 
the power law gives an infinite gradient at the wall. 
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Fig. la Meridional view of vane configuration. 

The primary convergence criterion was to examine 
the solution to see that the variables of interest were 
not changing with increased iterations. The primary 
variable of interest was surface heat transfer, because it 
converged at a slower rate than surface pressures. Be- 
cause of the large CPU requirements, most cases were 
not started from a simple initial flow field, but, when a 
parameter was changed, the previous solution was used 
as an initial guess. If the new parameter value pro- 
duced a small change in the flow field, the residuals did 
not necessarily increase significantly initially. With this 
approach, relying on “orders of maginitude” decreases 
in residuals was impractical. Low residuals relative to 
initial values without restart were verified, but were not 
the primary convergence measure. The inlet and exit 
mass flows differed by less than 1% at convergence. 


DISCUSSION of RESULTS 

Figure la shows a meridional view of the two configu- 
rations. They differ in the shape of the hub. For pre- 
sentation purposes, these configurations are labeled B 
for bellmouth, and S for the S-con toured hub. Config- 
uration B has lower acceleration near the front portion 
of the vane than configuration S for the same pressure 
ratio. The meridional view only indicates a part of the 
acceleration, since the flow is turned from social to an 
average of about 70 degrees. 

Figure lb shows the vane profile with the hub con- 
figuration B. The vane trailing edge angle varies from 
hub to tip. The highest amount of turning occurs near 
mid-span, and the lowest amount of turning at the hub 
endwall. This figure also shows the C-type grids used 
for the analysis. C-type grids were chosen in order to 
obtain good resolution of the flow field, and hence heat 
transfer, in the leading edge region. 



Fig. lb Typical grid used for analysis. 

Grid generation. For the test Reynolds numbers the de- 
sired near-wall spacing was 1 x 10" 5 of c x to obtain 
maximum values of yf of approximately one, which are 
needed to insure accurate results. Two approaches were 
taken for grid generation. A common grid was used for 
both flow analysis codes. Best results were obtained 
using an embedded grid with a non-matching condition 
along the cut line. A program, TCGRID, developed 
by Chima(1990), was initially used to obtain the grids 
used in the finite difference analysis. This program gen- 
erates a series of two-dimensional fully elliptical grids 
in the blade-to-blade plane based on Sorenson ’s( 1980) 
grid generation procedure. These grids, typically seven 
in number, are then stacked and interpolated to give 
a grid for each value of the spanwise index. For a 
177 x 49 x 65 grid, the seven blade-to-blade grids are 
interpolated to give 65 blade-to-blade plane grids. For 
a three-dimensional grid, the grid control parameters 
should be the same for each blade-to-blade section. Sat- 
isfactory three-dimensional grids with the desired near- 
wall spacing were not obtained using recommended or- 
thogonality constraint parameter values. Relaxing the 
parameters gave unrealistic flow solutions, probably due 
to excessive grid stretching. 

The approach finally taken for the finite difference 
analysis was to use the grid generation code developed 
for the finite volume flow analysis code, as described 
by Arnone et a!.(1992). In this approach the blade-to- 
blade grids are calculated in two steps. First a coarse 
grid, suitable for an Euler analysis, is generated using 
Sorenson’s(1980) technique. A grid with the desired 
near-wall spacing is then embedded within the coarse 
blade-to-blade grid. The resulting blade-to-blade grids 
are then stacked to yield the three-dimensional grid. 
Typically, the degree of orthogonality at the wall is less 
for an embedded grid, than for a fully elliptical grid. For 
turbine blade grids, the degree of orthogonality is signif- 
icantly improved if the constraint of matching grid lines 
along the cut line emanating from the trailing edge stag- 
nation point to the downstream boundary is removed. 
In contrast to the two-dimensional results presented by 
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Boyle and Ameri(1994), it was found that the conver- 
gence of the three-dimensional finite difference code, 
RVC3D, was significantly improved by not imposing a 
matching condition along the cut line. The finite vol- 
ume analysis, TRAF3D, of Arnone et al.(1992) was de- 
veloped to use grids which did not require a matching 
condition along the cut line. 

For each configuration the grid size was 177 x 49 x 
65. There were 56 increments on the pressure side of the 
cut line in the wake region. There were 32 increments 
along the pressure side of the blade, and 64 increments 
along the suction side. There were 24 increments on 
the suction side of the wake region cut line. 

Surface Pressures. Figure 2 shows a comparison of 
the predicted and experimental isentropic Mach num- 
bers, as determined from the static to total pressure ra- 
tio, for both configurations. Even though the passage 
geometry is complex, there is little spanwise variation in 
either the predicted or measured surface pressures. The 
computations are consistent with the experimental data 
in predicting the spanwise variation in pressures. These 
results were obtained using the TRAF3D code, and pre- 
dictions using the RVC3D code vrere nearly identical. 
Both CFD codes require the specification of either a 
uniform or average hub exit boundary pressure. Arts 
et al.(1990) have shown, for a similar vane, that when 
the exit Mach number exceeds unity, there was a con- 
siderable pitchwise variation in the midspan static pres- 
sure. At 0.42 chords behind the vane the maximum-to- 
minimum variation in static pressure was 11% of the 
dynamic pressure at a Mach number of 0.85, and in- 
creased to 23% at a Mach number of 1.1. The results 
shown in figure 2 are for a hub exit Mach number of 1.2. 
The blade surface pressure distribution was sensitive to 
the specified exit pressure. The tip exit pressure corre- 
sponded to an isentropic Mach number near unity, and 
had significantly less pitchwise variation than the hub 
exit pressure. Since only a few pressures were available 
to determine the appropriate exit boundary condition, 
the tip exit pressure was used. In the analysis, the hub 
exit pressure was adjusted so that the tip exit pressure 
matched the experimental value. 

Figures 3 and 4 compare predicted and experimen- 
tal hub and tip endwall isentropic Mach numbers for 
both configurations. Predictions were the same using 
either analysis, and the agreement with the data is 
good. Measurements, using approximately 53 taps per 
surface, extended from 0.06c x upstream of the lead- 
ing edge to 0.13c r downstream of the trailing edge. 
The predictions, which extend beyond the blade region, 
show the pitchwise variation in isentropic Mach number 
to be greater for the hub than for the tip. 




Fig. 2 Vane surface pressure distribution. 

Heat Transfer Comparisons. In this section pre- 
dicted and measured vane and endwall heat transfer 
are compared. Predictions of surface pressure distribu- 
tions were nearly identical using either flow analysis. 
However, predictions of surface heat transfer were dif- 
ferent between the two analyses. Since the analyses 
incorporated different turbulence models, these differ- 
ences could be due to either the analyses themselves, or 
in the turbulence models. It will be shown that the pri- 
mary difference in the heat transfer results between the 
two analyses is due to the choice of turbulence model. 

Comparisons of the predicted and measured vane 
surface Nusselt numbers are given in figure 5 using the 
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Fig, 3a Mis for configuration B hub. 

Baldwin- Lomax turbulence model, and in figure 6 using 
Chima’s turbulence model. Except for the suction sur- 
face transition region, there is little spanwise variation 
in the experimental Nusselt number. Even though the 
variation is not large, the results show that the midspan 
region of the suction and pressure surfaces have higher 
heat transfer than near either endwall. The peak suc- 
tion surface heat transfer occurs at transition, and is 
highest for the data at 50% and 98% of span. The lead- 
ing edge region heat transfer is relatively low compared 
to the average suction surface value. This results from 
the endwall convergence, which gives lower inlet veloc- 
ities compared with constant radius endwalls. 



Fig. 3b Mis for configuration B tip. 


The TRAF3D analysis agrees well with the data 
at the leading edge, s/c x = 0, for both configurations. 
However, this analysis overpredicts the pressure surface 
heat transfer prior to transition, which occurs at the lo- 
cation of minimum heat transfer. The RVC3D analysis 
underpredicts the midspan heat transfer at the leading 
edge, but is in good agreement with the pressure sur- 
face data prior to transition. The difference is a con- 
sequence of the TRAF3D prediction incorporating the 
Smith and Kuethe(1966) model to account for the ef- 
fects of freestream turbulence on laminar heat transfer, 
while the RVC3D prediction did not include a model 
for this effect. Without this model both analyses gave 
the same leading edge heat transfer rates. 
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Fig, 4a Afjs for configuration S hub. 

The predicted suction surface transition occurs 
downstream of the experimental location. This is a con- 
sequence of the way Mayle’s(1991) transition model was 
implemented. In this model the start of transition oc- 
curs when Res = 400Tu“ 5 / 8 . The measured upstream 
Tu was 6.5%. It was assumed that the fluctuations 
remained constant to determine the local Tu for the 
transition criterion. The upstream Mach number was 
0.14, so that the local calculated Tu at transition was 
approximately 1%. If the turbulence intensity is held 
constant, the predicted suction surface transition loca- 
tion occurs closer to the leading edge than is shown by 
the data. The results of Young et al.(1992), lead to the 
expectation that the primary effect of high freestream 



Fig. 4b Mis for configuration S tip. 

Tu was to cause early transition. They showed that 
the effects of high freestream Tu are at a minimum 
when either Res , or the length scale to boundary layer 
thickness ratio is large. Res is large on the endwalls, 
and the boundary layer thickness is small on the vane. 

Both predictions show the correct spanwise vari- 
ation in suction surface heat transfer after transition, 
where the highest values are at midspan and the low- 
est value are at 2% of the span. However, both analyses 
overpredict the amount of spanwise variation. This may 
not be due to defects in the turbulence model, since, as 
is shown in a subsequent figure, the vane suction sur- 
face heat transfer changes very rapidly near the endwall. 
The measurements may have averaged the heat transfer 
rates for a region close to the end wadi. 
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Fig. 5 Heat Transfer, Baldwin-Lomax(TRAF3D) 
model. 

The two predictions differ in the calculated heat 
transfer in the near wall region at 2 and 98% of span. 
This is especially evident near the leading edge for 
configuration S. The RVC3D prediction shows much 
higher heat transfer levels near the endwalls compared 
to midspan. The TRAF3D prediction shows smaller 
spanwise variation in this region. This is due to differ- 
ent assumptions made in the respective analyses. In the 
RVC3D code the eddy viscosity at a point is calculated 

from: fi t = yj(fit)\ 2 + (fit)? 2 - (fit ) i is found by consid- 
ering only the effect of the blade surface, and (fit ) 2 
found by considering only the effect of the endwali sur- 
face. Also, the length scale used for both values of fi t 
is the Buleev length scale, which can be approximated 
as the minimum distance to either surface. The 



Predctioo Exp. 



Fig. 6 Heat Transfer, Chima’s(RVC3D) model. 

TRAF3D code, however, weighs the individual contri- 
butions to the eddy viscosity. /i t = R(fit) 1 + (PtMl ” 
R). And, R = d\/(d\ + d\). For each value of fi t the 
length scale is the normal distance to the appropriate 
wall. In the corner region the eddy viscosity in the 
RVC3D code can be nearly 40% greater than the max- 
imum value from one surface. In this same region the 
TRAF3D assumption limits the value to the greater of 
the two values. Consequently, the RVC3D assumption 
gives greater eddy viscosity in the corner region, even 
for the same turbulence model assumptions. The re- 
sults in figures 5 and 6, show that the RVC3D code 
assumption regarding the eddy viscosity is in better 
agreement with the experimental data. 

Predictions were also made using the RVC3D code 
with the Cebeci-Smith turbulence model. Except for 
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the differences noted due to the implementation of the 
models, the Cebeci-Smith results were similar to the 
Baldwin- Lomax results shown in figure 5. 

Other differences in predicted heat transfer between 
the two turbulence models can be highlighted if transi- 
tion effects are absent. Figure 7 shows a contour plot of 
Nu as a function of span and surface distance for con- 
figuration B, and assuming fully turbulent flow. Simi- 
lar results were obtained for configuration S. This fig- 
ure shows comparisons of predictions for five different 
cases, and illustrates the effect of varying the turbu- 
lence model using the two flow solvers. The horizontal 
axis is the surface distance along the suction and pres- 
sure surfaces of the vane. Because of the endwall con- 
tour, the leading edge is at the location of maximum 
span. Part a of the figure shows Nu for the TRAF3D 
code with the Baldwin-Lomax turbulence model. Parts 
b, c, and d show RVC3D code results with Chima’s 
turbulence model, the Cebeci-Smith turbulence model, 
and the Baldwin Lomax turbulence model respectively. 
These three figures show the surface heat transfer for 
different turbulence model using the same flow solver 
and grid. The Cebeci-Smith model gives the highest 
vane heat transfer. The other two models give similar 
surface heat transfer. Comparing the heat transfer in 
parts a and d shows differences in heat transfer between 
the two flow solvers for the same turbulence model. 
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Fig. 7 Nu x 10 -3 for configuration B, fully tur- 
bulent assumption 

The differences are on the same order as the differences 
among the three turbulence models. The differences 
are not due to either differences in the implementation 
of the turbulence model, or to inherent differences be- 
tween the codes. If the fi t from the converged TRAF3D 
case was used in the RVC3D code, and the fit calcula- 
tion bypassed, the heat transfer was the same as in part 
d. The differences are due to the RVC3D code exhibit- 
ing a greater grid sensitivity than the TRAF3D code. 
Part e shows the vane heat transfer using a finer grid in 
the streamwise direction. Overall, this grid was about 
20% larger, but had a 50% increase in the number of 
grid lines, (from 32 to 48 increments), along the pres- 
sure surface. Sixteen increments were also added to the 
suction surface. The results in part e are in closer agree- 
ment with the part a results. When the TRAF3D anal- 
ysis was run with the finer grid, the vane heat transfer 
was neaxly identical to that shown in part a. 

In an effort to understand the differences among 
the results the manner in which the models were im- 
plemented was investigated. The TRAF3D code uses 
Sutherland’s Law to determine the variation of viscos- 
ity, and hence conductivity, with temperature. The 
RVC3D code uses a power law for this variation. Using 
the power law variation in the TRAF3D code did not 
change the predicted heat transfer. 

The Baldwin-Lomax model was applied everywhere 
in the flow field. In the RVC3D code a distinction was 
made, based on an input variable, between a “near wall” 
region and the entire flow field. For a j index greater 
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than a ./edge, the turbulent eddy viscosity was set 
to the outer region value, independent of where the 
crossover between the inner and outer regions occurred. 
Extending jedge to the grid size index resulted in no 
significant change in the predicted heat transfer^ Only 
when the turbulent eddy viscosity was suppressed com- 
pletely for j values greater than jedge was the heat 
transfer reduced by 10 to 15%. 

Both codes calculate the length scale in the same 
manner. For the vane or endwall the length scale 
is calculated as the normal distance from the grid 
line(RVC3D) or cell center(TRAF3D) to the appropri- 
ate surface. Taking the distance along the grid line for 
the length scale overestimates the length scale, and re- 
sulted in vane heat transfer rates up to 50% higher for 
both codes. This difference is due to the grids hav- 
ing significant degrees of non-orthogonality. Endwall 
heat transfer results were also sensitive to the approach 
taken to calculate the length scale. 

Figures 8 and 9 compare predictions of endwall heat 
transfer with experimental data for both analyses for 
configurations B and S respectively. The experimen- 
tal data extended from 0.11c r upstream of the leading 
edge to 0.11c r beyond the vane trailing edge. There 
were between 28 and 39 heat transfer measurements on 
each endwall. Upstream of the blade leading edge there 
was little variation in endwall heat transfer. The inlet 
Nu was calculated to be 1245. While the choice of the 
transition location noticeably affected the vane surface 
heat transfer, the location of transition on the vane had 
no noticeable affect on endwall heat transfer. Based on 
the measured inlet boundary layer thickness, the inlet 
endwall boundary layers was taken as fully turbulent. 

The predictions shown in these figures are for a 
177 x 49 x 65 grid. There was no significant difference 
in the predictions when a 209 x 49 x 65 grid, or a 177 x 
49 x 89 grid was used. 

Comparing figures 8a and 9a shows that configura- 
tion B had a somewhat lower hub heat transfer in the 
leading edge region, and a slightly higher heat transfer 
in the throat region. The predictions with any of the 
turbulence models do not show any significant differ- 
ences between the two configurations in the leading edge 
region. This is somewhat surprising, since, as is shown 
in figure 1, this is the region where there is a difference 
in the hub contour. In the throat region the Baldwin- 
Lorn ax model predicts the same hub heat transfer for ei- 
ther configuration, with a peak value that closely agrees 
with the measurements. Chima’s turbulence model un- 
derpredicts the peak heat transfer in the throat region, 
and is generally lower than the experimental data for 
most of the passage. The Cebeci-Smith model over- 




predicts the peak heat transfer in the throat region, but 
generally agrees well with the experiment ad data. 

In the tip region the experimental data are, as 
expected, relatively the same for both configurations. 
Again, there is a slightly higher throat region heat 
transfer for configuration B. The tip data show a region 
of high heat transfer at the top of the suction surface 
for configuration B, but not for configuration S. 
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4 





The predictions show the high heat transfer to be con- 
fined to only a small region of the tip endwall. Both the 
hub and tip predictions show the same relative agree- 
ment with the data for either hub configuration. 




The experimental data show a slightly higher tip end- 
wall heat transfer level compared to the hub endwall for 
both configurations. The predictions also show higher 
heat transfer on the tip endwall for a small region near 
the suction surface, close to the peak of the blade. This 
behavior is consistent with the experimental vane data 
shown in figure 5. At this vane surface region, 
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TRAF3D approach in the RVC3D code resulted in 
lower endwall heat transfer close to the suction surface. 

In an overall sense the Cebeci- Smith turbulence 
model tends to overpredict the endwall heat transfer, 
but to be in reasonably good agreement with the data. 
The Baldwin-Lomax results underpredict the heat 


transfer in the region near the suction surface, but 
agree reasonably well with the data elsewhere. Chima’s 
model underpredicts the heat transfer. This is espe- 
cially true for the maximum heat transfer in the throat 
region. 

Total Pressure Distribution. Figure 10 compares 
the total pressure distributions at 0.20 axial chords 
downstream of the trailing edge. Both computed re- 
sults overpredict the depth of the wake, 





Fig. 9 Nu x 10“ 3 for configuration S. 
but correctly predict the spanwise variation in the wake 
characteristics. The wake is thickest, and has moved 
further in the tangential direction, near midspan. Near 
midspan the minimum experimental total pressure ra- 
tio was 0.84. The RVC3D results were very similar 
for either turbulence model, and had a minimum total 
pressure ratio of 0.67. The TRAF3D results had a min- 
imum total pressure ratio 0.74. The spanwise variation 


Fig. 10 Total pressures 0.2c r behind vane, 
in the location of the minimum pressure shows that 
the parabolic variation in flow turning with span is 
well predicted. Also, both predictions are in agree- 
ment with the experimental data in that the wake is 
widest near midspan. The overprediction of the depth 
of the wake is consistent with the results shown by 
Boyle and Ameri(1994) in a two-dimensional analysis 
of the effects of grid topology on the aerodynamic and 
heat transfer characteristics of a high turning transonic 
vane. The tendency to overpredict the wake depth is 
not confined to the CFD methodologies used for this 
work. Results of the AS ME compressor rotor test case 
for Rotor 37(Strazisar,1994) showed that many analy- 
ses overpredicted the wake depth for the transonic com- 
pressor rotor. 
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PitehwtM averaged total preeaura ratio 



Fig. 11 Spanwise total pressure variation. 

The analyses predict a wake that is both too narrow 
and too deep. These effects compensate to some degree, 
and the pitchwise averaged total pressure downstream 
of the rotor is reasonably well predicted. This is il- 
lustrated in figure 11 where the spanwise variation in 
pitchwise averaged total pressure is shown. Both the 
experimental data and predictions were area averaged 
in the blade-to-blade direction to obtain these results. 
The experimental data extend from 6 to 94% of span, 
and show very low loss regions near the end walls. The 
maximum experimental loss occurs near midspan. Both 
predictions give too large a loss between the hub and 
midspan. The hub region is one of minimum static pres- 
sure, so that in the region between hub and midspan, 
overall loss is sensitive to shock losses. 

CONCLUDING REMARKS 

Heat transfer predictions were made using a finite 
difference(RVC3D) and a finite volume(TRAF3D) com- 
puter code for two vane geometries typical of actual en- 
gine applications. The test conditions of M 2 = 1.2, and 
Re 2 = 5.6 x 10 6 provided a severe test of the codes 
computational capability. Convergence was heavily in- 
fluence by the choice of grid topology. The grid gen- 
eration procedure described by Arnone et al.(1992) re- 
sulted in grids which gave satisfactory results for both 
codes. These grids were generated by embedding a 
near-wall grid. Primarily because of the required near 
wall spacing of 1 x 10“ 5 c x , it was found that using a 
non-matching condition along the cut line was highly 
desirable in terms of obtaining converged solutions. 

Solutions obtained with both codes converged in 
approximately the same number of CPU hours. Typi- 
cally, the grids used were 560,000 points, and the solu- 
tion was obtained in two hours on a Cray C90 computer. 


As expected the surface pressure distribution converged 
more rapidly than the surface heat transfer. Since the 
solution was stopped only after the heat transfer distri- 
butions were seen not to change over a large number of 
iterations, the actual convergence for either code may 
have taken somewhat less than two hours. 

Both computational analyses correctly predicted 
the spanwise variation in vane surface pressures, and 
the endwall static pressures for both experimental con- 
figurations. Best agreement with the data was achieved 
when the predicted exit static pressure was matched to 
the measured tip value. There was a larger pitchwise 
static pressure variation at the hub than at the tip. To 
minimize uncertainty, measured tip, rather than hub, 
values should be used for exit boundary pressures. 

Blade heat transfer predictions agreed reasonably 
well with the data for the pressure side of either vane for 
all turbulence model assumptions. On the suction sur- 
face of the blade Chima’s turbulence model produced 
results that agreed best with the data. This was inde- 
pendent of the choice of the transition model. The lo- 
cation of transition on the suction surface using Mayle’s 
model was bounded by the choice of local Tu used in the 
model. Assuming constant Tu resulted in early transi- 
tion, and assuming the fluctuation constant resulted in 
transition too far aft on the vane suction surface. 

Differences in predicted endwall heat transfer 
among the different turbulence models were greater 
than the differences in experimental heat transfer be- 
tween the two configurations. The results with the 
Cebeci-Smith turbulence model were in reasonably 
good agreement with the data, but tended to overpre- 
dict the peak heat transfer. The approach taken in 
the RVC3D code to determine in the corner region 
gave better agreement with the data compared to the 
approach taken in the TRAF3D code. There were sig- 
nificant heat transfer variations using the different tur- 
bulence models. Therefore, cases with significant varia- 
tions in Reynolds number should be examined to deter- 
mine if the degree of agreement with the data remains 
the same. This would increase the level of confidence 
in whichever turbulence model is used for predictions. 

Both computer analyses correctly predicted the 
spanwise variation in the wake behind the vane. Both 
analyses predicted too low of a value for the minimum 
total pressure in the wake region. This appeared to 
be due to underprediction of pitchwise mixing. The 
spanwise distribution of the pitchwise averaged total 
pressure was well predicted. The analyses and the data 
showed minimum total pressure loss close to both end- 
walls. 
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